attempt_single_step Interface

interface
public subroutine attempt_single_step(this, sys, h, x, y, f, yn, fn, yerr, k)

Arguments

Type IntentOptional Attributes Name
class(single_step_integrator), intent(inout) :: this

The single_step_integrator object.

class(ode_container), intent(inout) :: sys

The ode_container object containing the ODE's to integrate.

real(kind=real64), intent(in) :: h

The current step size.

real(kind=real64), intent(in) :: x

The current value of the independent variable.

real(kind=real64), intent(in), dimension(:) :: y

An N-element array containing the current solution at x.

real(kind=real64), intent(in), dimension(:) :: f

An N-element array containing the values of the derivatives at x.

real(kind=real64), intent(out), dimension(:) :: yn

An N-element array where this routine will write the next solution estimate at x + h.

real(kind=real64), intent(out), dimension(:) :: fn

An N-element array where this routine will write the next derivative estimate at x + h.

real(kind=real64), intent(out), dimension(:) :: yerr

An N-element array where this routine will write an estimate of the error in each equation.

real(kind=real64), intent(out), dimension(:,:) :: k

An N-by-NSTAGES matrix containing the derivatives at each stage.

Description

Attempts an integration step for a single-step integrator.